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ON A DISCRETE HILL'S STATISTICAL PROCESS BASED ON 

SUM-PRODUCT STATISTICS AND ITS FINITE-DIMENSIONAL 

ASYMPTOTIC THEORY. 
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» I ' Abstract. The following class of sum-product statistics 
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(where I, k = io and n are positive integers, < i < fc < n, P(p, h) is the set of 
all ordered parititions of p > into h positive integers and Yi^„ < ... < Yn^n 
are the order statistics based on a sequence of independent random variables 
Yl, y2,...with underlying distribution ¥{Y < y) = G{Y) = F(ey)), is intro- 
duced. For each p, T„(p)~^/f is an estimator of the index of a distribution 
j^ ' whose upper tail varies regularly at infinity. This family generalizes the so 

-4— J called Hill statistic and the Dekkers-Einmahl-De Haan one. We study the 

__j, limiting laws of the process {Tn{p), 1 < p < oo} and completely describe the 

covariance function of the Gaussian limiting process with the help of combina- 
torial techniques. Many results available for Hill's statistic regarding asymp- 
^ ^ totic normality and laws of the iterated logarithm are extended to each margin 

|/~v , Tn{p,k), for p fixed, and for any distribution function lying in the extremal 

QQ ' domain. In the process, we obtain special classes of numbers related to those 

VQ ' of paths joining the opposite coins within a parallelogram. 

p' 

cn 

o 

(N ■ 1. INTRODUCTION 

Let Xi, X2,... be a sequence of independent copies of a random 
X ■ variable (r.v.) with P (X < x) = F{x),x e M and F(l) = and let 

j^ I Xi„ < ... < X„„, denote the corresponding order statistics. Hill (1975) 

introduced 

(1.1) T„(l,fc,0 = ^ Yl j(l0g^n-i+l,n-l0gX„_,-„) 

(where /, k and n are integers such that 0< / = l{n) < k = k{n) < n, 
l/k — )■ 0, fc — )■ oo, k/n — )• 0, log stands for the Neperian logarithm) as 
an estimator for the exponent 7"^ > of a distribution function {d.f.) 
whose the upper tail varies regulary at infinity, that is : 

(1.2) VA > 0, lim (1 - F{\x) / (1 - F{x) = A"^ 

X— ^-|-oo 
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Mason in |T2j proved that the strong or weak convergence of T„(l, k, 1) 
for k = [n"] , < a < 1, to some 7~^ > for one value of a, < a < 1, 
is characteristic of rf./.'s statisfying (11. ip . This means in fact that 
F belongs to the extremal domain of attraction of a Frechet's law of 
exposent 7"^. It is well-known in extreme value theory that a d.f. F 
is attracted to some non-degenerate d.f.M, that is there exists two 
sequences of real numbers (a„ > 0)^>q and (fen)„>o such that. 

(1.3) Vx G M, lim F"(a„x + 6„) = M{x), 

iff M is the Frechef type of d.f. of parameter 7 > 0, 

(T A\ ( \ — j exp(— x"^) X > 

^ ' ^ '^^ ^ 1 elsewhere ' 

or the Weibull type of d.f. of parameter 7 > 0, 

(1 5) ih (x) = I ^^P(~(~^)'^) a; < 

^ ' ■' r-y\ ) \ -^ elsewhere 

or the Gumbel type of d.f., 

(1.6) A(x) = exp(-e"^), x eR 

A huge amount of statistical applications for modeling extremes ares 
based on these limiting laws. For instance, Resnick (1987) cited a wide 
range of phenomena described by extreme value theory. Further, in 
order to characterize the d.f.'s attracted to one of these three d.f.'s, L6 
(see |T^) introduced 

(1.7) 

^ i=k 

Tn{2,k) = T X] ^(l-'^ii/2)(logX„_i+i,„-logX„_j,„)(logX„_j+i,„-logX„_j,„) 
1=1+1 

and showed that the couple (T„(l, k, /), T„(2, k, I)) combined with some 
other auxilliary statistics characterizes the whole domain of attraction 
and each of the three domains of attraction. These characterizing sta- 
tistics are studied in details in [8] and [9], corresponding laws of the 
iterated logarithm are given. Tn{2,k,l) may be viewed as a second 
dorder form of T„(l, fc, /). In the same view, Dekkers, Einmahl and De 
Haan introduced a moment estimator in [5] 



'1-8) "^" ^ 7: X^ z(logX„_i+i,„ - logX, 



n—i,nj 
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And it happens that T„(2, k, 1) = 2A„. The importance of this second 

order estamator is that the couple (T„(l, k, /), T„(2, k, /)) separates the 

whole domain of attaction in the sens that the couple (T„(1, k, I), (T„(1, k, /)T„(2, k, /)~^^^)) 

takes three limiting values according to three the possibilities : Fg 

D{(f)^), Fg D{ip^), Fg -D(A). Theses previous works of many authors 

cleary suggested extensions to higher orders and advocated a stochastic 

process view. Our aim is to generalize results available for T„(1,A;,/) 

and T„(2, k, I) for higher order forms T„(p, k),p > 1, and to study the 

limiting behavior of the process {T„(p, k),l < p < +00} . 

Let us make further notations before defining Tn{p,k,l) = Tn{k). 
Let P{p, h) be the set of all ordered partitions of p > into positive 
integers, 1 < h < p : 

(L9) 

P{p,h) = {{si...Sh),yi, I < i < h, Si > 0; si + ... + Sh = p} ,1 < h < p. 

Let Yi^n < ••• < Y'n,n be the order statistics of the n first of a se- 
quence of independent random variables Yi,Y2,... such that P{Y < 
y) = G{y) = F(e^), y >0. Now define for 1 < / < /c < n, p > 1, io = ^, 
(LIO) 

P to ih-1 






h=l(si Sh)&P(p,h) -11=1+1 ih=l+l 



It is easy to check that T„(l) is Hill's statistic, T„(2) is that of (1.3) 
and, for instance, for p=3, we have 



-. j=k . i=k i=j . 

J-ny^j) = T\ / ^ 7 {yn~j+l,n ~ ^n-j) + / ^ / ^ 7^l[^n-j+l,n ~ J^n-j,n) 
j=l+l j=l+l i=l+l 

^ \^ n—i+l,n ^n~i,n) ' \^ n—j+l,n ^n—j,n \^ n—i+l,n ^n—i,n) 
i=j m=i 

~r / ^ / ^ ^^ \J-n—m+l,n ^ n—m,n) \^ n~i+l,n ^n—i,n) ^ 
i=Z+l m=Z+l 

(Lll) 

l.-'^n— j + l,n ^ n—j,n) i 



GANE SAME LO. 



Our aim in this paper is to establish the asymptotic normal- 
ity theory for the finite-dimensional distribution of this process. Our 
best achievement is the convergence of the finite-dimensional distribu- 
tions of the process T„ = (T„(p), 1 < p < oo) when suitably centered 
and normalized to those of a Gaussian process {/(p), 1 < p < c>o} . The 
computation the covariance function of this Gaussian process required 
much combinatorial calculations. This combinatorial work, along with 
the study of the numbers classes which appear for the varaince calcu- 
lation, is the one of the main part of this paper since the results to 
be used for the d./.'s F are largely developped in [S], [H] and [lU] and 
other authors. 

First, we state limit theorems in section II along with the methods to 
compute the covariance function of I{p)- Section 3 is devoted to prove 
the combinatorial relations in section 2. Finally, section 4 is devoted 
to the proofs or the theorems. The reader is is referred to [7], [H], pT] , 
[13], and [6] for general references on extreme values theory. 



Before we state our results, we mention that a continuous general- 
ization of the Hill statistics has been given and studied ([5]). 

2. STATEMENT AND DESCRIPTION OF THE RESULTS 

We begin with the description of a class of numbers involved in this 
work. We introduced three classes of numbers that will be used for the 
covariance computations. 

Definition 1. The non-negative integers (3{v,r), v > 0, r > 1 defined 
by 

i. Vf > 0, /3(f , 1) = 1 

ii. Vf > I, /3{v,2) = 1 
iii. Vr > 3, P{1, r) = (3{2, r - 1) + (1, r - 1) 
iv. Vr > 2, /3(0,r) = /3(l,r-l) 

V. Vt; > 2, r > 3, (3{v, r) = (3{v + l,r - 1) + f3{v - 1, r) 

are called the type I numbers. 



We have the two following clog rules associated with points (iii) 
and (v) of this definition. 



V 


u-l-v 






V 


u 






u+v 






u 
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The repeated application of them leads to the computation of all the 
numbers of this class. For instance for v < 10, r < 10, we have the 
different values of /3(z/, r) 



v\r 




2 


3 


4 


5 


6 


7 


8 


9 


10 







1 


1 


2 


5 


14 


42 


132 


429 


1430 


1 




v = 1 


u + v = 2 


5 


14 


42 


132 


429 


1430 


3862 


2 




u = 1 


3 


9 


28 


90 


297 


1001 


2432 


7294 


3 






4 


t = 14 


48 


165 


572 


2002 


6072 




4 






5 


t + s = 20 


75 


275 


1001 


3640 






5 






s = 6 


27 


110 


429 


1638 








6 






7 


35 


154 


637 










7 






8 


44 


208 












8 






9 


54 














9 






10 
















10 






1 

















Remark 1. The second clog rule is that of the numbers f3*{v, r) of paths 

in 1? joining (0, 0) to {y + r, r) within the parallelogram [(0, 0), (f , 0), (f + r, r), (r, r)] 

which is 



(2.1) I3{v,r) 



v + 2r ;-M l^l \,r<v + l 



(See p^). The differences between /3*(-'-) and /3(., .) are the following: 

(1) (f,r) satisfies r < v -\- 1 for /?*(.,.) while f > 0,7 > 2 are 
arbitrary in /?(.,.). 

(2) The first lines {v = and v = 1) differ for /3 * (., .) and /3(., .). 

Definition 2. The non-negative intergers iir{0,v,6), I < S < t verify- 
ing 

i V(5> 0, Vr > 1, /i^(0,t;,r) = 1 
ii V(5, 1 <(5<r, ;u^(0,0,(5) = 1 

iiiV(5, l<(5r, Vt; > 1, /i^(0,t;, 5) = /i^(0, w - 1, 5) + (0,t;, 5 + 1) 
are called the r — class of type II numbers. 



Remark 2. Point (Hi) of this definition yields the clog rule 



u 




u+v 


V 



which enables to compute easily these numbers . Here are some exam- 
ples. 
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r 



v\6 






v\6 


1 


2 




^;\5 


1 


2 


3 










1 




1 


1 


1 




1 




1 


2 




2 


3 


2 




2 




2 


3 




3 


6 


3 




3 




3 


4 




4 


10 


4 




4 




4 


5 




5 


15 


5 




5 




5 


6 




5 


21 


6 





Finally, 

Definition 3. The r — class of type III numbers are the non-negative 
integers fir{^,v,6), v > 0, t > 1 satisfying points (i), (Hi), (iv), (iv) 
anf (v) of definitionUl for 6 = r and, at the place of (ii). 



v+l 



(2.2) 



V(t;>l), /i,(l,i;,2) = 5^/i,(0,A;,l) 



fc=i 



Here are examples of computation 

r = 1 



r 



v\5 




2 


3 


4 


5 




v\8 




2 


3 


4 


5 







1 


2 


5 


14 







1 


5 


14 


42 


2 




3 


9 


28 




1 




5 


14 


42 


132 


3 




4 


14 






2 




9 


28 


28 


90 


4 




5 








3 




14 


48 






5 












4 




20 




















5 













Tab 2.5 



Tab 2.6 



We finish by remarking that these tables are very quickly filled with 
standard softwares. We are now able to describe the gaussian process 
involved here. 

2.1. THE LIMITING GAUSSIAN PROCESSES. 

Definition 4. The time series {I{r),r = 1,2, ...} is called an extremal 
Gaussian process if and only i/E(/(r)) = for r > 1 and its variance 
and covariance E(/(r)^) = (T^(r) and E(/(r)/(p)) = a{r, p) statisfy for 
some positive functions Ci(.) and C2(., ) 



i) a(0) = 1, a(l) = 2 



3=r 



ii) V(r > 1), a'^ir) = Ci{r-)a{r) and a{r) =2 E/5(l,i) a{r - j) 
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j=r 

iii) VI < r < p, a{r,p) = C2(r) ^/ip_^(ll, j)(a(r - j) 

i=o 

with by convention ^,-(1,1,1) = /iT-(0, 1, l),/iT-(l, 1, 0) = 1 for all 
r > 1. 

Here again, the values of o"^(r) and cr{r,p) are easily computed 
forCi(r) =Ci{r,p) = 1. 



4 


5 


11 


29 


70 


3 


4 


9 


20 




2 


3 


6 






1 


2 








r/p 


1 


2 


3 





Tab 2.8 



This table cleary shows that this process is not not stationary since, 
for instance. 



(2.3) 



a(2, 1) = 3 7^ 9 = a{3, 2) ^^ 29 = a(4, 3) 



Definition 5. The time series {I{p) + e{p)Z, p = 1,2,...} ,where e > 
1 is a real function and Z is a standard Gaussian r.v. such that E,[I[p)Z) = 
— 1 for all p = 1,2, ..., is called a reduced extremal process.. 

2.2. LIMIT THEOREMS FOR T„(p).. Let |/o = sup {a;, G{x) < 1} , 
0< X < z < yo and 



(2.4) 

and for p> 2 

(2.5) mp{x,z) 



mi{x, z) 



;i - G{t)) dt 



(1 — G{t)) dt dyi...dyp_i 

'X Jyi -Jyp-i 

with mp{x,yQ) = mp{x), p > 1. In the remainder of the paper, we 
shall use, without any loss of generality, the following representations 
of the order statistics of Y : 
(2.6) 

{n,n < ... < yn,n,n> 1} = {G-\Ui,n) < G-\U2,n) < ■■■■ < G-\Un,n),n> 1} 



GANE SAME LO. 



where f/i,„ < ••• < Un,n are the order statistics of a sequence of inde- 
pendent uniform r.w's on (0, 1). Put now 

(2.7) 

x„ = G-\l-k/n),x^ = G-\l-Uk+i,n), Zn = G-\l-l/n),Zn = G-i(l-f/z+i,„) 

and 

71 
(2.8) rp{Xn,Zn) = -mp{Xn,Zn),P>l 



Finally, let D{(pj),'y > 0, be the set of (i./.'s F satisfying II. 21 D{ipj) 
the d.f.'s such that Xq = sup{x,F{x) < 1} < +oo and F{xo — ^) E 
D{(p-y) and -D(A) the set of d./.'s F such that 

j.^ F-\l - xu) - F-\l - u) ^ logx 
^■^ «->o F-i(l-yn)-F-i(l -m) logy 

for all X > 0,y > 0,y ^ 1. 

It is clear that T = D{A) U D{if) U D{ip) where £>(</?) = U^>oD{(f^) 
and -D(V') = U^>oD{4'^), is the set of all d.f.'s attracted to some non 
degenerate d.f. We have the first theorem limit. 

Theorem 1. Let F G F, I be fixed and k satisfy 

(K) < k = k{n) — )■ +00, k{n)/n — > an n — > +cxd, 

then {A;-*^/^ (^n(p) — 7p(a;„)) /rp{xn), I < p < +c>o} , converges in distri- 
bution to the extremal gaussian process process in the canonical topol- 
ogy of N°° with 



(2.10) Ci(r) = n i ^^^ }> , r > 1 

and 



(2.11) C2ir,p)=n\ l^l \,r>l,p>l,r^p, 

for < 7 < +00. 

Before giving the next theorem, recall that Fe F may be represented 
by constants c and d and by functions f{u) and b{u), < u < 1, with 
b{u) and f{u) tending to zero as u tends to zero, through 

(2.12) G~^(1-m) =logc-(logM)/7+ / b{t)t-^dt, < m < 1, 

J u 
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For F G D{<^^), 

(2.13) yo-G-\i-u) = c{l + f{u))u^/'^exp{ I b{t)r^dt),0 < u < 1 

J u 

and for F G D{;ip^),a.'n.d 

(2.14) G'\l -u) = d- s{u) + / s{t)r^dt), < m < 1, 

J u 

for F G D{A), where s{u) = c(l + /(«)) exp(/J 6(t)rMt), < m < 1. 
(I2.12P and (I2.13P are the Karamata representations while 02.14p is the 
de Haan-Mason-Deheuvels one. 

Replacing Tp{xn) by the non-randon sequence rp(x„) requires regu- 
larity conditions on f(.). Such characterizations for p = 1,2 are given 
in L6 (1991a). In fact, they will hold again. But since we are only 
interested in putting reduced processes to the fore, we can only use the 
simplest condition, that is : f ' has a derivative in some neighborhood 
of zero and 

(RC) lim ufuQ = 0. 

u — >0 

We are now able to formulate our second theorem. 
Theorem 2. Let F eT, I be fixed, k satisfy {K) and {RC) hold, then 
{/i;^/^(r„(p) - Tp{xn))/rp{xn), l<p < +00} 

converges in distribution to the reduced form of the extremal process 
of Theorem U\ in N°° induced with its canonical topology with e{p) = 

(7 + P)/7; < 7 < +CX). 



Remark 3. Following the notation in L6 (1990b), the case < 7 < 
+00 corresponds to F E D{ip^), 7 = +00 means F G D{A)UD{(f). Un- 
less the contrary is specified, the corresponding values for the functions 
in 7 are obtained, in the second case, by letting 7 — > +00 in the first 
case. For example, for F G D{K)UD{ip), e{p) = hm^^+oo(7+p)/7 = 1- 



Remark 4. These two theorems prove the existence of the gaussian 
processes introduced in Definitions\J\ and\E and show how they may be 
observed and simulated. 
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3. Technical lemmas concerning numbers generation 
Recall that for all j > 1, mj{x,yo) < +oo for F G F, 

h]{l) = dpi dpi / dp2 / ... / '"^ ^^, mj{qr-i,Zn) dqr-i, 

Jpo Jqa Jql J qi -^Ppr-2 

where f G N, j G N*, r > 2, qo = pi, po = Xn ', 



v\ 



(3.1) 

-f]{0,6)= dqi dp2... ''"'^"^^ , "^^ mj{qr-s+i,Zn) dqrr-5+h 

J x„ Jqi J Qt-S ^' 

and 

(3.2) 



dqs-1 
' qi •JqT-i J Xn J qT J v\ •^Ps-2 



rz„ pzn pzn Mt rz„ pzn rq 

7]'(1,5) = / dqi dp2... / dqr / dpi / dqr+i / ... / 

J x„ J ai J aT-^ J Xn J q-T J v^ J vx 



(3.3) / < r 5 1 U_ (^^^^_^^^^)l^^^^^_^, 

where z/ G N, r G N*, (5 = 2,3,... 

The three class of special numbers given here appear when com- 
puting these integrals. 

3.1. COMPUTATION OF h^^(r). 

Lemma 1. For all j > 1, r > 2, z/ > 0, the ratios 

(3.4) I3iiy,r,j) = h''j{r)/mj+^+2{r-i)ixn, z^) 

are positive integers and depend only on {u, r) so that 



(3.5) Vz/>0, V>2, /3(z/,r,j) = /3(z/,r). 
Proof. First put r = 2. We have 

rZn PZn /„ r) )'^ 

(3.6) /i^"(2) = / dpi ^^' /'' m,{qi,Zn)dqi, v > l,j > 1. 



PI ^! 



By remarking that. 



(3.7) mj{qi,Zn) = -dmj+i{qi,Zn) /dqi, 

and by integrating by parts, we arrive at 
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(3.8) h'<{2) = h^^-li2), for u > 1, j > 1. 
But, 

(3.9) /i°(2) = mj+2{xn, Zn), for j > 1, 
so that 

(3.10) /l5(2) / mj+o+2{2-l){Xn, Zn) = P{0, 2) = 1. 

By repeating (13. 8p until its right member becomes h^{2) and by using 
dil, we get 

(3.11) V(j > 2),V(z/ > l),h'<i2) = m,+,+2{xn,Zn) 
which proves the statements of the lemma for r = 2, that is 

(3.12) V(j > l),V(z/ > 1),1 = /i,^(2)/m,+,+2(2-i)(a;„,^„) = /3(z/,2) 

Now, by applying again a change of variables like (13. 7p and by inte- 
grating by parts in /ij(7), we obtain 

(3.13) V(j > 1), V(r > 3), Vz/ > 1, /i^^(r) = h'^Xl{r - 1) + h^-\{r) 

By assuming that the statements of the lemma hold for r — 1 > 2, that 
is 

(3.14) V(j > l),V(z/ > 1), /3(z/,r-l) = /i,"(r-l)/m,+,+2(r-2)(a:„, z„), 

and by repeating (I3.13P until the second term of its right member be- 
comes an /i2(7) term, we show that the expression /ij(r)/mj+i,_|_(r-i)(x„,, Zn) 
/3(z/, 7) does not depend on j for any z/ > 1 and that 

h=u-l 

(3.15) V(r > l),(Vz/>l), /3(z/,r-l) = ^ /3(z/ - /i,r - 1). 

h=-l 

This proves that /3(z/, r) is also integer for z/ > and t > 3. This 
together with (13.121) proves lemma by induction. It is easy to derive 
from this latter that 

(3.16) V(r>3),V(z/>2), /3(z/,r) = /3(z/ + l,r - 1) + /3(z/ - l,r). 
Further, for r > 3, J > 1, we obviously have 

(3.17) /,0(^) = /,i^^(,_i) 
so that 
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(3.18) V(r > 3), Vj > 1, /3(0, r) = /3(1, r - 1) 

By applying (13. 7p and by integrating by parts, we arrive at 

(3.19) V(r > 3), h]{r) = h]_,,{{r - 1) + h]_,,{^), 
which, combined with (I3.17P and (I3.18p . imphes 

(3.20) V(r > 3), /3(1, r) = /3(2, r - 1) + /3(1, r - 1) 
Finally, set by convention 

(3.21) V(z/> 0), /3(z/,l) = 1 

and see that fl3J6D . fl3:T8|) . fl3:20|) and fl3:2T]) together show that the 
beta nimbers are generated by the ratios 

(3.22) K<{r)/mj+y+2{T-i){xn,Zn), for j > 1, z/ > 0,r > 1. 

D 
3.2. COMPUTATION OF YA'^,^)- 

Lemma 2. For all 5, 1 < 5 < r, z/ > 0, the ratios 



; ^ ^ ^5 



(3.23) 7j(0, r)/mj+^+^_5+i(a;„, 2;„) = /i^(0, z/, 5,j) 

are integers depending only on r, z/ and 5 so that 

(3.24) 

Vr > l,V(5,0<(5<r,Vz/>0,Vj > 1, /i^(0, z/, 5, j) = /i^(0,z/,5) 

Proof. A (13.70 - like change of variables yields 

(3.25) ^-(0,r)=7;;ii(0,r), j > 1. 
It is easily checked that 

(3.26) 7°(0, (5) = mj+r-S+l{Xn, Zn), 

which gives 

(3.27) W S,l<5<T, ^r{0,0,5) = l 
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Combining this with f l3.25p . we get 

(3.28) Vz/>0,/i^(0,z/,r) = 1 

By using also a (I3.7p -hke change of variables, we have 



(3.29) 

V J > l,V(z/ > 0), V(l < (5 < r), 7^(0,(5) = 7-^^(0,5 + 1) +7)7i'(0, '5) 

One proves this lemma by induction over r—r, r = 0, 1, ..., r— 1 through 
(I3.29P after having taken into account the initial column given in (I3.28p . 
This induction yields 

(3.30) Vz/> 0,V1 <6 < T,fir{0,'^,S) = ^,/i^(0,z/,r), 

fc=0 

which, in turn, gives the following clog-rule 



(3.31) Vz/> 1,V1 <5<r, flr{0,U,S) = flr{0,U-l,6)+fir{0,U,S + l) 

It is now proved that the r — class of numbers are generated by the 
integrals 7^(0,(5). D 

3.3. COMPUTATION OF 7j(0,(5),5 > 2.. 

Lemma 3. Let t eW be fixed. For all j > 1, v > and 5 > 2, the 
ratios 

(3.32) 'yj{l,6)/mj+^+r+2{S-l)iXn,Zn) = fir {1, 1^,5 J). 

are integers depending only on (z/, 6) so that 

(3.33) 

\freN*, V(j>l), V(z/>0), V(5>2), fi,{l,u,6,j)=fir{l,u,6) 

Proof. :First, we compute 75^(1, 2). A (I3.7|) -like change of variables gives 

(3.34) 7j(l,2)=7Si'(0,l)+7;7i'(l'2), J>1,^>0 
By repeating this latter u times, and by using 

(3.35) Vj>l, 7?(l,2)=7l+i(0,l), 
and, finally by applying Lemma 121 we get 
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k=iy+l 



(3.36) 7j(l, 2) = ( ^ /i^ (0, k, 1) mj+u+T+2 {xn, Zn) ■ 

fe=l 

This proves that the statements of the lemma hold for 5 = 2 with 

(3.37) Vz/ > 0, /i^ (1, z/, 2) = ^/i^ (0, fc, 1) 

fc=i 
From now, this proof follows the lines of that of Lemma [T] with exactly 
the same methods. We also get the same conclusions. D 

4. PROOFS OF THEOREMS. 

Csorgo-Csorgo— Horvath and Mason have constructed, in [IJ, a proba- 
bility space {Q,U, P) carrying a sequence of Brownian bridges {Bn{s), < s < 1} 
n = 1,2.... and a sequence of independent uniform r.v.'s on (0,1) 
Ui, f/2, ... such that for all < z/ < 1/4, 



(4.1) 


\^(U^(s) - s) - B^is)\ 

i<s<l-i (S(l-S))i-Z/ 

n — — 71 ^ 


= Opin-" 


and 






(4.2) 


\^{S-Vn{s))-Bn{s)\ 
l<s<l-i (S(l-S))i-Z/ 

n — — n ^ 


= Op(n-": 


where 







J 
Un{s) = - for Uj^n < S < f/j+l,n, < S < 1, 

is the uniform empirical distribution function and 

Uj,n for '-^<s<i-,l<j<n, 0<s<l 



^"^'^ " ^ K(0) = f/i, 

is the uniform quantile function and finally, and 

are the order statistic of t/i, ...f/„. 

Throughout these proofs, we suppose that we are on this proba- 
bility space. Consequently, the sequence Yi,Y2,... defined above and 
the sequence of empirical distribution function based on them, will be 
represented as 
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(4.3) 

{Yi,n, l<i<n,n>l}= {G-\l - f/„-.+i,„), l<i<n,n>l} 

and 

(4.4) 

{1 — Gn{x), — oo <a;<a; + oo,n> 1} = {Un{l — G{x)), — oo < x < +00, n > 1} 

First, routine calculations yield 

(4.5) Tn{p) = - / ... / f/„ (1 - Git)) dt dy,...dyp^^,p > 1 

^ J Xn Jyi J Vp-l 

The details of the computations are omitted. The reader may verify 
it for p = 1, 2 and 3. Now, let 

(4.6) an{s) = n^'\Un{s)-s),Q<s<l- 



(4.7) mi,p(x„)=/ / ... / Bn{l-G{t))dt dyi....dyp_,; 

Jxn Jyn Jyp-1 

pXn rXn pXn 

(4.8) m2,p{zn)= / ... / Bn{l-G{t))dt dyi....dyp.i- 

^ Zn J Vx 'J Vp-l 

•2n f^n 



ms,p{xn,Zn)= / ... / Bn{l-G{t))dt dyi....dyp-i] 

J Xn Jyi J Vp-l 

(4.9) 

.^s 1 rXn r-Xn rXn 

^"=(^) / / - / (^n{l-G{t))-Bn{l-G{t)) dt dyi....dyp.^; 
and 

(4.10) 

Wn{p) = (ly r r ... r B,, (l - G{t)) dt dy,....dyp^. 



'Xn Jyi Jyp-1 

We have 



Lemma 4. For all p > 1, 

1 

(4.11) Vk{Tn{p) - Tp{Xn,Zn)) = W„(p) + -R„ + f t] ^ "^l,p(Jn) 
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"n\2 sr^ [Zn-XnY /- N 1 Sr^ mi^j[Xn) , ~. 

X I T ) > . -; m2,p-j(2n) + 7T > . -; m3,p-j (X„, 2n) • 



Proof. This is straighforward. D 



We already have all the necessary row matierals (to handle the error 
terms of Lemma |4] ) in Lemmas 4.1 and 4.2 in [9], summerized as 
follows. 



Lemma 5. let F ^T with F{1) = and et G be associated with F by 
G{x) =F(e^),x > 1. Then 

1- Rp{x,F) ~ (xo — xY {^''jZ^iil + j)^^} , as X ^ Xq, whenever 
F G Di^/j,) 

2- Rp{x, G) ~ i?i(x, Gy, as x^ yo, whenever F e {K)UD{^) 

3- {{z — xY / Rp{x, G) — )■ +00 and Rp{x, z, G) / Rp{x, G) — )■ 1, as 
X -)■ xo, -2 -)■ xo, (1 - G{z)) / (1 - G{x) -)■ 0. 

Much details on how using Lemma [5] to treat errors terms in Lemma 
m are given in the proofs [9] . We therefore omit them only for sake of 
conciseness. We finally get. 

Lemma 6. Let F G F. Then for all p > 1, 



(4.12) 

Vk{Tn{p) - Tp{Xn)) I TpiXn) = W'n(p) / Tp^Xn) + Op{l), aS U -^ +CX) 



Also, we have 
Lemma 7. Let F G F and (RG) hold. Then, 



(4.13) 

^/k{Tp{Xn) - Tp{Xn)) / Tp{Xn) = e{p) uk'^^^ iUk,n ) + 0^(1^ 



-e(p) i?„ ( - ) + Op(l), as n -)■ +cx) 
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Proof. (Outline of the proof). 
Check that 



(4.14) Tp{Xn) - Tp(Xn) = -{Xn - X„)mp_i(x„) + /?„(2) 

n 
where 



-|-Rn(2)| < \Xn -X„|^SUp ( -,1 - G{Xn) ] +^ " ., " TUp^ji 

n \n J j^ r- 

By Lemma 3.6 in [8] and Lemma [5] below, 



Xr, 



T + 1 A; 

(4.15) Vk(Xn - Xn) I Rl(Xn) = BJ-) + Cp(l) aS Tl -^ +00 

7 n 

It is not difficult to show that \/kRn{2) / Tp(xn) -^p as n — > +oo 
by using Lemma |5] and (I4.15P and that f (1 — G{xn)) — t- 1 as n — t- +cxd 

. Thus, Vk(Tn{p)-Tp(Xn,Zn)) / Tp^Xn) and \/k{Tn{p)-Tp{Xn)) / Tp{x) 

behave asymptotically as 



(4.16) Wn{v) / rp{xn), P>1, 

and as 



(4.17) Wn{p) I TpiXn) - e{v)Bn{-) 

n 

U 



Let us say a few words on the finite-dimensional distributions (f.d.d) 
of (I4.16P and (I4.17p . Since W„(p) is a multiple Riemannian integral 
for each p > 1, it is clear that, for n fixed, any linear combination of 
iWnipi), ..., Wn-{pj)) , for any fixed j > 1, is limit everywhere of linear 
combinations of f.d.'s of the Brownian bridge Bn{»-)- But (14.160 and 
(I4.17P are well-defined normal random variables and then, their f.d.d 
are gaussian. We have now to verify that their limiting covariance 
functions are finite. 
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4.1. COMPUTATION OF THE VARIANCE OF Wnir),r > 1.. 
Recall that h(s,t)=E (5„(1 - G{t) (1 - G{s)) = min 1-G{t), (1 - G{t))- 
(1 - G{t)) (1 - G{s)) , for < s, t < 1. One has 

(4.18) V(r > 1), E {Wn{r)) = 

and 



(4.19) 



n2^ ^ 



n / -^n / ^n / ■^n / ^-■ri 



Vr>l, E(H^„(r)^) = - / //•••/ / h{s,t) ds dt 
xdpi..dpr-i dqi... dqr-i 






(4.20) =-/ / H{p,,q,)dp,dq,= :a^(r). 



n '^ 't-n 



Since for all {p,q) {p < yo, 1 < Vo), H{p,q) = H{q,p), cutting the 
integration space into (pi < gi) and (pi < gi) yieds 



(4.21) al{r) = 2- dp J H{p,,q^)dp, 

1^ Jxn Jxn 



One has (see [2] and [9]) 
(4.22) 



n / -^n 



yixn<y<Zn), / / Ms,t)dsrft = 2 / / l-G(s) rft 1 (l+r„)(l) 
J?/ Jj/ ^Jj/ Jt 

with |r„(l)| < 1 — (^(xn). For the remainder, we shall proceed by 
induction. Suppose that for r > 2, 

(4.23) 

V j, (1 < j < r-l),Vx„ <y < Zn, / (igi / dpi / (ig2 / dq2...dpj-i / rfs. 

/i(s,t) rfs rft = 2a{i) m2jiy,Zn) (1 + r„(j)), 

where |r„(j)| < 1 — G{xn) and a(j) does not depend on the size n. 
Now, by cutting 

■ dpj+i 
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into 

■dpj+i + / -dpj+i 
'Pj Jii 

in fl4.2ip sequentially for j = 1, ...,r — 1 and by applying (14.231) . we 
get 



(4.24) al{r) = 2^ Y.^{r - j) h\^_^^^, (j) (1 + r^{j, r)) 

where |r„(j, r)! < 1 — G{xn) and, by convention, a(0) = 1 and h\{l] 

1 X m2+=l{Xn,Zn) 



Now, by applying Lemma [H 



3=r 

(4.25) al{r) = (2^ ^a(r - j)/9(l, j) (1 + ^nW) m2r{xn,Zn) 
with |r„(r)| < 1 — G(x„). This and fl4.22p prove that for all r > 1, 

(4.26) al{r) = 2- a{r) m2r-(a;„, 2;„) 



Now, according to Lemma [5], we get 



j=r 



(4.27) hm E(iyO {{ry / r^ixn)) = 2Ci(r) Va(r - j)/3(l, j) 

where Ci(7) = 1 for F G D{A)UD{ip), 

Ci(7) = n (^^1 /orF G D(^,). 

4.2. COMPUTATION OF THE COVARIANCE FUNCTION 

OF Wn{-)- One has, r < p, r = p — r, 



ly-. r ^n r Zn P^n /* Zn f* Zn 



(4.28) 
E(H^„(r) Wn{p)) = 't I dqi I ... / dq^ I I dpi I dqr+i 



n I ^n 



dq2 / dqr+2 / / h{s,t) ds dt 
Cutting f^" •Avj+i into , f^^ •.dpi + f^" •.rfpi and P" .dp,- into 
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, /* .dpj into J^^+' • dpj + J^""^^ • dpj, we get 

(4.29) E{W^{r)W^{p)) = - {a{r - l)72V-i)+i(0, 1)(1 + r^{r, 1)) x 

r 

^(a(r-j)72V-i)+i(l'.^')(l + '"n(r,j)) 
i=2 

where |r„(r, j)| < 1 — G(x„), for 1 < j < r. By Lemma |3l 



(4.30) 



3=r 



E{Wn{-f)Wn{p)) = < ^^^^('"-j) /^r(l,l,j) > (l + r„(r,p) m^+p(x„, z„) 

with |r„(r, p)| < 1 — G{xn) and by convention, /^^(l, 1, j) = /^r(0, 1, j). 
Now by Lemma 



(4.31) Mm E{Wnir)Wnip)) / r,(x„) rp(x„)) = C2(r,p) a(r,p), 

n— >+oo 

where for 1 < r < p, 

j=r 

(4.32) a(r,p) = ^a{r - j) pp_^(l,l,j) 
and /or F G D{ip^), 

(4.33) C2(7,p)=n^ ^ + ^ 



i=i 1 7 + P + J 
and for F e D{A) UD{^), 

4.3. COMPUTATION OF THE COVARIANCE FUNCTION 

OF {wn{r)/Tr{x„)} — e{r)Bn{-). It is immediate that 

(4.34) 

E{Wn{r) Bn{-)) = {1 + rn)Ymr{xn,Zn), where |r„| < 1-G{xn). 
n k 

Put Wn{r) = {lV„(r)/r„(r)} - e(r)S„(J). Thus, 

(4.35) hm EiWnir) W{p)) = C^ir, p)a{r, p) - e(r) - e(p) + e{r)e{p) 



A DISCRETE GENERALIZATION OF HILL'S STATISTIC. 21 



while 



(4.36) lim EWniry = 2Ci{r) a{r) - 2 e(r) + e(r)^ 

n— >oo 

We have now finished to compute the covariance functions oiWn{r)/Tr{xr. 
and Wn{r) via (lOTj) . ( lOTj) . (lOSj) . and (gSS]). These formulas say 

that ?i'n(^)/Tr(a;n) is an extremal gaussian process and Wn{r) is the 
reduced process. 

These processes deserve a great interest. Particularly, the structure 
of the covariance function should be investigated. High order proper- 
ties, as p — )■ +00, must be obtained. All that will be done further. 

We just mention now that laws of the iterated logarithm in [9| may 
be easily extended to T„(p) as follows : 

(4.37) 

-^ n \P) ^p \ •^n ) ^n j 

Vp > 1, lim sup(resp. inf .) — - — = p{p) ( resp.—p{p) ), 

n^oo Tp{Xn)v2k loglogU 



(4.38) 



-'-nXP) 'TpV^n) cls 



Vp > 1, lim sup(resj9. inf .) — ■'' = a{p) ( resp.—a{p)), 

n^oo Tp{Xn)\/2k loglogU 

where a{p) = a'^ip) — e(p)(2 — e(p)). 
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